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Abstract. Using the renormalisation group (RG) we study two dimensional 
electromagnetic coulomb gas and extended Sine-Gordon theories invariant under 
the modular group SL(2,X). The flow diagram is established from the scaling 
equations, and we derive the critical behaviour at the various transition points of 
the diagram. Following proposal for a SL(2, Z) duality between different quantum 
Hall fluids, we discuss the analogy between this flow and the global quantum Hall 
phase diagram. 



1. Introduction 

In statistical physics, self-duality in the sense of Kramers and Wannicr maps the 
high temperature regime of a given model with its low temperature regime. In two 
dimensions, examples of self-dual theories are provided by the 2 (Ising),3 and 4 states 
Potts, the Ashkin Teller, and the clock models. All of them can be represented in 
terms of an electromagnetic Coulomb gas (Kadanoff 1978, Kadanoff 1979), allowing to 
understand this self-duality as an exchange of the electric and magnetic components of 
the charges of the equivalent Coulomb gas (Kadanoff 1978). This self-duality allows 
to locate the transition point, whose study requires however a full renormalization 
group study, cither directly on the Coulomb gas formulation (Nienhuis 1987) or on 
the associated Sine-Gordon theory (Boyanovsky 1989). 

More than fifteen years ago, Cardy uncovered a generalization of this simple 
duality by adding a topological coupling between magnetic and electric charges of a two 
flavors Coulomb gas (Cardy 1982). Motivated by the study of oblique confinement in 
four dimensional Z p lattice gauge theory with a topological 9 term, Cardy & Rabinovici 
(1982) formulated a Coulomb gas where the presence of the 9 term considerably enlarge 
the usual Kramers Wannier duality g — > l/g to the full modular group SL(2, Z) (Cardy 
1982). This 9 coupling was later extended to higher dimensions in (Shapere & Wilczek 
1989). By analogy with the work of Kramers and Wannier, Cardy (1982) derived the 
location of the numerous transition points in the presence of this topological coupling, 
as the invariant points of the modular group. 

The purpose of this article is to explicitly study using the renormalization group 
the scaling behaviour of an extension of the Cardy-Rabinovici model. Besides the 
critical behaviour associated with the transition points identified by Cardy, this 
renormalization study allows to find the scaling flow of the model : the fixed points 
defining different phases are defined and related to the transitions points. We thus find 
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a condition that two different phases must fullfill to be related by a transition. This 
is, to our knowledge, the first explicit renormalization of a two dimensional modular 
invariant model. 

Besides the applications of this extended Coulomb gas to the previously cited 
statistical models in two dimensions and to fermion models in 1+1 dimension via 
its Sine-Gordon representation, such two dimensional modular invariant theories are 
of interest for the study of the global phase diagram of the quantum Hall effect. 
This effect correspond to the quantization of the transverse conductivity a xy of a two 
dimensional gas of electrons in a strong transverse magnetic field, and it is now well 
understood in terms of the microscopic Laughlin waves functions (Laughlin 1983). The 
problem of the nature of the transitions between the different quantum Hall fluids is of 
particular interest and remain unsolved. In particular the notion of superuniversality 
was proposed to explain the similar behaviour of all the transitions in the global phase 
diagram of the quantum Hall effect (Kivelson et al. 1992). Such a superuniversality 
can be deduced from a duality of the underlying model, relating all the transition 
points with each other. Indeed the similarity between the phase diagram of (Cardy 
1982) and the expected renormalization flow diagram in the two parameters scaling 
model of the Quantum Hall effect was noticed (Shapere & Wilczek 1989) before this 
proposal. Later on Liitken & Ross (1992) related more precisely the properties of this 
phase diagram with the presence of a SL(2, Z) symmetry, in the framework of a two 
parameters a xx , a xy scaling theory. 

On the other hand, the experiments of (Shahar et al. 1996) on the transition 
between the v = ^ and the Hall insulator have shown a reflection relation between 
the nonlinear current density and electric field on both side of the transition. Within 
the effective description of the quantum Hall effect , this reflection was interpreted 
as a duality which exchanges the electric and magnetic component of the low lying 
excitations (Shahar et al. 1996). Motivated by these experimental results, several 
authors focused on modular invariant models (Fradkin & Kivelson 1996, Pryadko & 
Zhang 1996). Even though no derivation of a microscopic modular invariant model 
exists, these studies focused on the general constraint from the modular symmetry on 
the phase diagram and expressions for the conductivity at the transitions (Liitken & 
Ross 1992, Dolan 1998). The phase diagram of the modular invariant Coulomb gas 
we study in this paper is expected to mimic the one of the quantum Hall effect. As 
an explicit renormalization of a modular invariant model in two dimension, we expect 
this study to be hclpfull for the description of the critical behaviour of the quantum 
Hall transitions. 

In another context, it is interesting to notice that modular invariant Coulomb 
gas also appeared in the work of Callan & Freed (1992) on the dissipative motion of 
a charged particle in two dimension in a transversed magnetic field and a periodic 
electric potential. Using a mapping to a one dimensional version of the model of 
(Cardy & Rabinovici 1982), these authors found a phase diagram as a function of 
the dissipation strength and the magnetic field which resembles the diagram obtained 
from our renormalization study. 

This article is organized as follows : in the section 2.1 we define the model both 
in its Sine-Gordon version and as a lattice Villain model with topological term. 
The mapping to an electromagnetic Coulomb gas is then obtained (2.2) together 
with the symmetry SX(2,Z) (2.3). In the next section (3.1) we derive the general 
renormalization group equations for the model using a Kosterlitz scheme, and the 
phase diagram is established in part 3.3 together with the critical behaviour at the 
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critical points of this diagram. 



2. Electromagnetic Coulomb gas 

2.1. Extended Sine- Gordon model and theta terms 

In this article we will consider the two dimensional extended Sine-Gordon model 
defined in its more general form by the partition function 

a •> 

with the action 

2k 



/A2 / 1 
E [g-Td^d^ 

\ a, 6—1 

+2y cos(0 a ) + 2yJ2 ™sp (b" 1 ]^ 6 - M«V ) ) (1) 

a a / 

In this definition we consider 2k integer valued fields (j) a and their dual fields </>° 
defined by idi<f> — eijdj(j), where i = 1, 2 labels the two directions of the plane, and dj 
is the antisymmetric tensor. The coupling matrix g ab and M ab are 2/c x 2k matrices; 
M ab has to be antisymmetric for renormalisability of the model (see the following) : 
M ab = -M ba . In this context the operators 

C?„, m (r)=cxp i (^n a ^(r)+^m a a (r)j (2) 

\ a a / 

correspond to the parafermions operators of Fradkin & Kadanoff (1980) which create 
a vector electric charge n and magnetic charge m in site r. 

Alternatively we can consider a Z p gauge theory defined on the square lattice by 
2k Villain models perturbed by symmetry breaking fields of strength y and coupled 
by a (topological) term : 

AviUain Iff" 1 ] " ~ 2 ^<<) - 2? TP<i) + M») E « 

a 2—1,2 a. 

E <c - ^ E Mafc ^- - 2n P A h) - 2n p<j) ( 3 ) 

a a 

where the A% i correspond to the integer valued gauge field defined on the bonds of 
the square lattice, and d x <i> a a is the discrete derivative 4>% +x ~ H ere and in the 
following summation over repeated indices is assumed. 

The meaning of the coupling between the different fields <fi a becomes clearer upon 
restriction to a two components model {k — 1) with the coupling constants 

g ac = g6 ac . M ac = ^_ £ ac ^ = y ^ 

2tt 

By considering the 4 component fields (j)^ — (0, 0, c/) 1 ,^ 2 ) which depends only on x\, X2 
we can write the quadratic term in (3) as F^F^ U where v = 1, ... 4 and is the 
usual electromagnetic tensor associated with the field ^ : F^ — d^ v — d u (f>^ — 
2npA^ u . the second term can be interpreted as a coupling between electric particle 



Renormalization of two dimensionnal modular invariant Coulomb gas 



4 



and the field <^ while the last term reads {0 /2-K)e^ pa F^ FP a = (e/2n)F^F^ v , which 
is known as a topological 9 term. This model and its Coulomb gas formulation in two 
dimensions was first studied by Cardy & Rabinovici (1982). The general M coupling 
in (1) is thus the analog of this topological coupling for the 2k real components field 
in two dimensions. Other extensions to higher dimensions of space were developped 
in (Shapere & Wilczek 1989). 

Although in the following, we will derive renormalization group equation for the 
general model (1), we will only analyze in details the scaling behaviour of the simpler 
model (4) with only electric charges in the first component (a = 1) and magnetic in 
the second (see below). This model is defined by the following restriction of (1) : 

H= 4^g~J 5 : (W 1 ) 2 + W 2 ) 2 )+y cos ^ 1 )+y cos p(.^ 1 ^ 2 -^^ 1 ) (5) 

We can now express the partition function of the above models in terms of 
electromagnetic Coulomb gas, extending the usual case of (Nienhuis 1987). 

2.2. SL(2,Z) invariant electromagnetic Coulomb gas 

To obtain the Coulomb gas formulation of the above model, we first use the Villain 
approximation (Villain 1975) of the cosine coupling in (1) : 

e 2ycos0 ^ e in.<f>+n 2 ln(y) 

n=0,±l 

This approximation is valid for small coupling strength y. The models with both 
forms of the interaction are known to be in the same universality class without the 
topological term (Boyanovsky 1989), which ensures the same in our case. Within this 
approximation, the partition sum of (1) consist of a trace over the fields a (r) and the 
electric and magnetic charge density n a (r) and m a (r) of the exponential of the action 

y a.b a a 

+ ln( ? y)EK) 2 +lnyE(™ a ) 2 ) 

a a / 

In this bare model the vector charges n(r) have component 0, ±1, however upon coarse 
graining charges with higher components will be generated. The lattice Villain model 
with a topological coupling (3) leads exactly to the same sum, where the magnetic 
charges m a (r) are located on the sites of the dual lattice. They are defined by 
the oriented sum of the potential A a over the plaquettc surrounding the dual site 
: m a (r) = -e^A^. 

After integration over the field (f> a (r) , and using the neutrality of the charges 
/ n a (r) = J r m a (r) = (imposed by the infrared regularisation) , we obtain an 
electromagnetic two dimensional coulomb gas of electric n a and magnetic m a charges, 
defined either on a lattice (and its dual for the m charge) or in the continuum, which 
take value in Z 2fc . It is defined by the grand canonical partition function 

Z= E ( fi / J e^M (6) 
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where the y njm are the charge fugacities, and the primed sum counts each distinct 
configuration only once. The corresponding action can be written in its most general 
form as 

A cg [n, m] = i [g ac (n a + P M ab m h ) * G * (n c + P M cd m d ) + p 2 (g- 1 ) ac m a * G * m c ] 

-i pn a * $ * m b (7) 

where we assumed summation over repeated indices and used the contraction notation 
: n * G * n = J2i^j n iG{ri — Vj)nj. In this action, the potentials G and $ are defined 
respectively by the propagators 

9 ab G(r) = (rm h (o)) - (r(0)cf> b (r)) 

-i$(r) - (</> a (O)0 b (O)) - (<f> a (0)4> b (v)) 

In a neutral Coulomb gas, these propagators need only to be regularised at short 
distances : in the following we will use a real space hard cut-off a, corresponding to 
hard core charges. The asymptotics of these interactions is given by G(r)+«$(r) ^ r >a 
ln(x + iy)/a. For a definition of these potentials on the lattice, see (Kadanoff 1979, 
appendix A). 

2. 3. Symmetries for two components charges (k = 1 ) 

We now considere more precisely the model (5). With the presence of the 9 (or 
M) coupling, the usual Kramers- Wannier duality g <-> p 2 g~ 1 of the electromagnetic 
coulomb gas (Kadanoff 1978) is considerably enlarged (Cardy 1982). Besides the time 
reversal symmetry T S[{n a }, {m a }, g, 9] = S[— {n a } 7 {m a }, g, — 6], the 2tt periodicity 
of the 9 coupling translates into P S[{n a }, {m a }, g, 9} = S[{n a - e ab .m b }, {m 6 }, g, 9 + 
2tt]. Finally the action (7) is also invariant under the self-dual transformation 
D S[{n a },{m a }, 5 ,6>] = S[{e ab .m b },{e ab .n b },g\9'} with 

| = ^W ; g0 = - 9 '0' (8) 

These transformations are better parametrised using the complex coupling constant 
z = p-§^+ipg~ 1 ■ the above transformation simply read P(z) = z+l and-D(z) = — 
As obvious from this writing, D and P do not commute, and they generate the whole 
infinite discrete group SL{2, Z). The purpose of this letter is thus to study explicitly 
the behaviour of the modular invariant Coulomb gas (7) under the renormalization 
group. 

3. Renormalization a la Kosterlitz 

3.1. Renormalization group equations 

Without a 9 (M) term, the electromagnetic Coulomb gas (7) can be renormalized either 
by following the (Anderson- Yuval) Kosterlitz scheme (Nienhuis 1987) or directly by 
an operator product expansion in the Sine-Gordon formulation (Boyanovsky 1989). 
Both methods amounts to study the product of the parafermion operators (2) and give 
the same results. Hence for seek of simplicity, we will follow the Kosterlitz approach, 
extending the classical method to the presence of the 9 term (see (Nienhuis 1987) for 
a detailed review of this method). 

Upon coarse graining the model by increasing the hard-core cut-off a, we leave 
the partition function (6) invariant by defining scale dependent coupling constants 
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and fugacities. Three different contributions to these variables have to be considered : 
naive rescaling, fusion and annihilation of electromagnetic charges. The naive rescaling 
conies simply from the change of cut-off in the integration measure and the interaction 
G(r) : its gives the eigenvalue of the fugacities y n , m - Upon infinitesimal increase of 
the cut-off a — > a = ae dl , the distance between two neighbors charges can become less 
than the new cut-off a. When these two charges form a dipole, we integrate them out 
(annihilation of charges), while we simply glue them into a single charge otherwise 
(fusion of charges). In this last case we get a correction to the fugacity of the new 
charge, which together with the naive rescaling, reads the following scaling equation 
for the fugacities : 

%n, m = ( [ 2- g -Y,{n a +pM ab .m h f - |! ^( m «)^ ^ m ( 9a ) 

+71" J]] <5n'.m"+n".m'2/n',m'2/n",m" 

(n' .m^ + fn" ,m") — (n,m) 

with the notation n.m = ^2 a n a m a . 

In the other case, the small dipoles we annihilate upon coarse graining screen the 
interaction between distant charges. In the partition function (6), the term involving 
two (non zero) charges (n s , m t ) and (n s , m t ) distant from a < \r s — r t | < ae dl can be 
expanded in a/|r, — | (the distance between distant charges) and yields 

-2ir 2 dl ^2 ya i , mi yn j , mj {Q.jsQ.ks ~ Pjs/3ks)yl 3 ,m 3 G ( r i - r fc) 

J^i-r^a ( n . >m .) ; ( n .. ;m ..) 

where a js = g ac NjN^ + p 2 (g~ 1 ) ac m'*m c s , fij s = p(n°-m a s + n£m") and we have defined 
the composite electric component of the charges N? = n a - + p M ab m b . Using 
the antisymmetry of M ab and reexponcntiating this contribution, we can check the 
renormalisability of the model to order y 2 as this contribution can be cast into 
a contribution to the matrices g, M and the fugacities. We obtain the following 
corrections to g and M to the order y 2 : 

d l9 ac = -27T 2 ]T (g ab g cd N b N? - p 2 m a s m c s )y 2 ^ ms (9b) 
= -2n 2 p ]T [(g-^Kmt (g^rNtmt] yl s ^ e (9c) 

n s ,m s 

Notice the antisymmetry of the correction (9 c) to M ab . We can now derive the 
renormalization flow from these scaling equations. 

3.2. Specific model and charge asymmetry 

In the following we will restrict our study to the model (5). Contrarily to the usual 
case where the fugacities of electromagnetic charges are symmetric both in n and m 
(Nienhuis 1987), here this symmetry is broken by the presence of the 8 coupling. As 
seen above with the T symmetry, changing n to — n amounts to change also 6 to —8, 
and similarly with the transformation m — > —to. The only symmetry which does not 
modify either g or 9 corresponds to (n, m) — > (— n, —to), which ensures the neutrality 
of the gas for any value of the couplings. Hence we cannot assume the to or n parity 
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of the fugacities y n , m as the condition y n . m — y- n , m = Vn-m wm n °t be preserved 
by the RG. 

For the model (5) we need only to consider charges satisfying n 2 = 0,-m 1 = 0. 
We will use the notation (n, m) for (n 1 ,?™ 2 ). The above RG equations (9) can then 
be written in a more pleasant way using the notation x = pg^ 1 an d t = p 8/2n : 

8 lX = 2tt 2 P ]T [(n + tmf - x 2 m 2 } yl m + 0(y 3 ) (10a) 

n,m 

diy n , m = 2 (l - £ (n + tmf - |zm 2 ) y n , m + 0(y 2 ) (10c) 

In the last equation the non linear term from (9 c) has been forgotten, being 
subdominant around each transition point. However it has to be taken into account 
when analyzing the topology of the whole phase diagram. The above RG equations 
correspond to the starting point of our analysis of the different transitions. Their 
invariance under the duality transformations is made more explicit if we write them 
as a single equation for the complex coupling constant z : 

diz = i2ir 2 pJ2( n + zm ) 2 yl,m (11) 

n.m 

On this last equation, the invariance under duality of the flow (z — > — z^ 1 ) is obvious. 

The invariance under duality of these beta functions, which was assumed in 
previous studies of modular invariant models, have important consequences : as the 
axis t — is invariant (in the limit of vanishing fugacities) and thus correspond to a 
flow line of the above equations, we know from the action of the modular group on 
this flow line that all the flow lines of the phase diagram are either straigth lines or 
arcs of circles. The flow diagram obtained by numerically integrating these equations 
is shown on figure 1, where only charges with m < 3 have been taken into account. 



3.3. Analysis of the flow 

To derive the phase diagram from the renormalization analysis, we first find the 
domains where a charge (n, m) proliferate. Such a charge proliferate, or equivalcntly 
the parafermion operator (2) associated with a charge (n + j^m, m ) is relevant, when 
the associated fugacity increases under rescaling. From the renormalization eigenvalue 
of the fugacities in (10c) we deduce that the pure electric charges (n, 0) are relevant 
only for x > pn 2 /4. For x < 1 the composites charges (n, m) become relevant inside 
the circles [n,m] centered on (t,x) = (—n/m,2/(pm 2 )) of radius 2/{pm 2 ). All this 
circles are thus tangent to the axis x = (see figure 1). 

Each circle [n,m] is caracterised by a ratio v = n/m instead of n, m, as all charges 
(An, Am) with A e Z are generated under renormalization in this circle (see figure 
1). In the following we will write v — n/m where m is the minimal magnetic charge 
allowed in the circle. Under renormalization, in a circle v — n/m, x = pg^ 1 flows to 
zero while t is renormalized to t* = v. Hence each rational point of the axis x = 
correspond to a fixed point of the RG, caracterising a given phase. The phase m = 0, 
stable for x > pn 2 /4, can be viewed as a special circle of infinite radius, caracterised 
by v = oo. In this phase x — > oo and t is slightly renormalized to a (non universal) real 
value. Below x < pn 2 /A and between all these circles, all the fugacities renormalizc 
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Figure 1. RG flow in the plane p^,pg^ x . The flow are shown up to m = 3. 
The regions inside the circle corresponds to stability regions for charges (n, m) 
with nonzero n, m labelling the size of the circle while n gives their positions. For 
pg^ 1 > 1 only charges with m = exists. The transition points between m ^ 
phases arc shown as black points while the black region corresponds to a phase 
with no charge (neutral phase). 



to zero, which correspond to a neutral phase of the Coulomb gas, caracterised by non 
universal renormalized g and 9 and vanishing fugacities. 

Transitions in this phase diagram arc a priori of two different types : the 
transitions between two circles caracterised by different ratio v and the transition 
between a phase with relevant charge (n, m) and the neutral phase. Transitions of the 
first kind correspond to tangent points between circles in the phase diagram. Such 
points exist only when p = 4. For p < 4, the circles overlap and these transitions are 
no more accessible by the present perturbative study, while for p > 4 all circles with 
different ratio v are disconnected and these transitions disappear. In the following we 
will consider in more detail the case p = 4, whose phase diagram is shown on figure 1. 

A transition between two phases (ratios) v\ and v 2 is allowed if the two 
corresponding circles are tangent to each other, which can be written 

Vl-V2 = ±—^~ (12) 
miTO2 

The corresponding transition point is located in (t, x) = {— (n\m\ + n 2 m 2 )/(m\ + 
m|),I/ \m\ + 777*2)). To analyze further these transitions, we need to remark that 
any circle v can be deduced from the circle v = by succesive application of the 
transformations P, P^ 1 and the duality D. Moreover the circle v = itself is the 
image under D of the line of phase transitions x = 1 . The transition points between 
the phase v = and other v' = ±l/m (see figure I) are all images under D of the 
transitions points between [1,0] and the phases [77, 1], which themselves are related to 
the transitions [1,0] <-> [0, 1] by action of P and P~ x . Thus all the transition points 
between two ratios v and v' correspond to the same critical behaviour. Similarly 
all transitions between a phase v and the neutral phase (no charge condensate) are 
images of the transition between [1,0] and this phase. We can now study these two 
transitions perturbatively in the y n ^ m . 

First at a transition between u\ = 771 /mi and v 2 = n 2 /m 2 we can write 
renormalization equations for the distance orthonormal to the transition point Si 2 
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while the distance perpendicular to this axis is not rcnormalized to lowest order. 
Using e = (m\ + m%)8 12 and y = y nimi , y = y n2m2 we obtain 

die = 2n 2 (y 2 - y 2 ) d t y = 2ey d t y = 2ey 

which correspond to a Kosterlitz-Thouless like diverging correlation length £ ~ 
exp(cte/|e| 2 ). In both phases, the correlation functions of the parafermion operators 
decays exponentially : 

e i(n0(r)+m0(r)) e -i(n0(o)+m0(o))\ ^ g -| 



The transition between the neutral phase and a circle v is describe by a critical 
point with one marginal direction (tangent to the circle) and the same diverging 
correlation length. While the correlation function have the same kind of exponential 
decays from the circle side of the transition, they decay algebraically in the neutral 
phase : 



e i(n0(r)+m0(r)) e -i(n0(o)+m0(o)) 



4. Analogy with the Quantum hall effect 

Within the context of the two parameters scaling theory of the quantum hall effect, 
the first idea consists in identifying the coupling of the Coulomb gas model with 
the components of the conductivity tensor of the electron gas a xy = t, 
The ratio v defined in the above study naturally translate into the filling factors 
of the quantum Hall states : each circle of the flow diagram correspond to a given 
quantum Hall fluid and the renormalization equations we obtain provide the expected 
asymptotic value for this conductivity : a* y — v, a* x — 0. Following this analogy we 
can find that transitions between two filling factors v\ , v 2 are allowed if they satisfy 
the expected relation (12). All the allowed transitions between two plateaux are in the 
same universality, being related by the duality D (8) which exchange the electric and 
magnetic charges, as in (Shahar et al. 1996), thus supporting the superuniversality 
hypothesis (Kivelson et al. 1992). 

However we notice that in our flow diagram, filling factors with any denominators 
exist, while only even denominators are expected for quantum hall fluids. Moreover 
fixed points corresponding to filling factors (circles) with even and odd denominators 
are related by modular transformations : hence one cannot avoid the even 
denominators in this model. This result has to be related with recent works (Georgelin 
et al. 1997) which reveal that the symmetry hidden behind the quantum Hall 
hierarchies should be the subgroup T(2) of SX(2,Z) instead of the whole modular 
group. This subgroup is generated by the transformation z — » — z/(2z+ 1) (instead 
of z — > — z^ 1 in this study) which preserves the oddness of filling factors. This 
transformation, when restricted to the axis t = (imaginary axis), correspond to 
g — > g/(l + 4g 2 ). Hence it does not exchange the high and low temperature phases 
of the usual electromagnetic Coulomb gas(Kadanoff 1978). Thus an extension of the 
usual electromagnetic model to a T(2) invariant theory will not simply correspond to 
the addition of a topological term as for the model studied in this article. 

To conclude, let us comment on the comparison with the recent studies which 
focus on the constraint on the beta function from the modular symmetry (Dolan 
1998). This study (together with (Burgess & Lutken 1997)) assume a modular (or 
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T(2)) symmetry in a model with two parameters a xx , a xy , and derive the general form 
of the beta functions of this model. In our work the model (6) possess an infinite 
number of coupling constants corresponding to the electromagnetic charge fugacities. 
Thus a direct comparison between the beta functions is not possible. However it should 
be noticed that our renormalization procedure provide the first explicit example of the 
commutation of the modular symmetry and the renormalization, which is the central 
hypothesis of (Burgess & Lutken 1997, Dolan 1998). 

5. Conclusion 

In this article we thus extended the renormalization study of (Nicnhuis 1987) to the 
Coulomb gas with a topological 9 term in two dimensions. The renormalization scheme 
we used provides scaling equations which were found to be themselves invariant under 
the modular group. To our knowledge this study provide the first example in two 
dimensions of an explicit renormalization of a modular invariant model. The rich 
phase diagram of this model certainly deserves more work : of particular interest 
will be the applications of this renormalization study to the extensions of the various 
statistical models, related to the electromagnetic Coulomb gas in the absence of the 
9 term. 

Note added : after completion of this work, a preprint (G. Cristofano & Nicodemi 
1998) appeared on a related quantum Hall model where a finite electromagnetic 
background charge is added to the k — 1 restriction of the model (7) . Although 
these authors used a simplified scheme, their scaling analysis seem to agree with the 
complete RG equations we derived in this article. 
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